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Quantum information is an exciting field promising a revolution in information 
processing. A key ingredient for the advancement of this field is the development 
of technologies that can implement quantum information processing (QIP) tasks. 
This paper describes recent progress using nuclear magnetic resonance (NMR) as 
the platform. The basic ideas of NMR quantum information processing are detailed, 
examining the successes and limitations of liquid and solid state experiments. Finally, 
a future direction for implementing quantum processors is suggested, utilizing both 
nuclear and electron spin degrees of freedom. 

PACS numbers: 

Quantum information processing is changing our fundamental understanding of what 
information is and how it can be manipulated. Recent work has lead to experimental 
proof-of-principle demonstrations of the control of small quantum systems, and new devices 
capable of harnessing larger quantum systems could change the technological landscape of 
the 21st century. There exist many proposals to physically realize such a quantum processor 
and presently a few of these models are able to manipulate quantum bits (qubits) : quantum 
systems whose observables are given by the Pauli matrices. A qubit encodes a fundamental 
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unit of quantum information. One such proposal uses the nuclear and electron magnetic 
moments and is the subject of this short review. 

Many atoms nuclei possess a magnetic moment. When placed in an external magnetic 
field, these moments result in discrete energy level systems that can be manipulated with 
resonant electromagnetic radiation leading to NMR. It was first observed more than half 
a century ago by Purcell and Bloch [HE]- It has become a powerful analytical tool with 
many applications such as non-destructively determining molecular structures for chemistry 
as well as static and dynamic imaging in both industry and medicine. 

Many nuclei such as or ^■^C are spin-1/2 quantum system. These are ideal qubits. 
The initial state of the system is obtained by allowing it to thermalize. By irradiating these 
nuclei at the appropriate frequency it is possible to rotate their individual states one at a time 
leading to the so called single-qubit gates. Nuclei affect each other through interactions that 
can also be controlled, leading to two-qubit gates. By composing these two sets of gates we 
can reach any unitary transformation. This is known as universality [3]. After an algorithm 
consisting of a series of these gates, the final state can be observed by measuring the current 
induced by the rotating magnetic moments of the sample in a conducting coil. Algorithms 
are to be performed on a timescale shorter than the characteristic decoherence time of the 
system, which is about 100 times the gate time in liquid state NMR. 

The strength of NMR technology is the exquisite control that can be implemented in 
multi-qubit systems. Many decades of radio-frequency (RF) engineering improvements today 
provide sufficient control to implement experimental benchmarks such as quantum error 
correcting codes and simulations of quantum physics systems, among many others. In this 
paper, we first give a brief review of both liquid and solid state NMR and then discuss recent 
work using both electron and nuclear spins to form more powerful information processing 
devices. The latter work provides a pathway to reach high purity states which has been a 
weakness of the liquid state NMR proposal. 



The rest of the paper is organized as follows: sections [TAj- |IE| give a detailed intro- 
duction to the concepts of NMR based quantum information processing in the liquid state 



environment; sections |I F| - |I G| introduce more advanced ideas of quantum control based on 
refocusing and optimal control methods; section |TT] describes the solid state environment, a 
three-qubit solid state processor and experimental results demonstrating high-fidelity con- 



trol; section III discusses a new direction involving systems of coupled electron and nuclear 



spins, including motivations, prospects, and recent progress; and finally section [TV] concludes 
the paper. 



I. LIQUID STATE NMR 

A. Magnetic interactions 

In the semi-classical picture, the spin of a nucleus behaves like the dipolar moment of 
a magnet possessing angular momentum parallel to its magnetic moment. When placed 
in a constant magnetic field pointing along a certain direction, (customarily defined as 
the z direction) the dipolar moment precesses around this axis. The frequency of this 
precession is called the Larmor frequency and is dependent on the external magnetic field, the 
nuclear isotope and its chemical environment within the molecule. For quantum information 
purposes, we are mainly interested in spin-1/2 nuclei (e.g. ^H, ^^C, ^^N, ^^F, ^^Si and ^^P 
to name a few). 

Placed in magnetic fields generated by modern superconducting magnets, different species 
of nuclei have differences in Larmor frequency on the order of MHz. For example, the Larmor 
frequency of is about 500 MHz in a 11.7 Tesla magnet, while that of ^^C is about 125 
MHz. Depending on the symmetry of the molecule, two nuclei of the same species can either 
have the same Larmor frequency, or can have a frequency difference (called chemical shift) 
ranging from a few Hz to several kHz. Typical liquid state NMR experiments involve an 
ensemble of around 10^° identical molecules dissolved in a solvent whose effect on the nuclear 
magnetic moments of our molecules can be neglected. 

When two spins are spatially close, their dipolar moments interact with each other. The 
strength of this coupling is dependent on the distance between the two spins and their relative 
orientation with respect to the external magnetic field. In a liquid, the molecules move and 
rotate around each other on a much shorter time scale than the interactions occurring 
between them. This causes the intermolecular and intra-molecular dipolar interactions to 
average to zero on the NMR time scale (i.e. the Larmor period time scale). In solid state 
NMR, however, dipolar interactions remain and can be controlled as discussed in section |Tl} 

Within the same molecule, there are still interactions between the spins in the liquid 
state. If the wavefunctions of bonding electrons overlap spatially with a pair of nuclear spins, 
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the electron mediates an effective interaction between the nuclear spins. This interaction 
is independent of the external magnetic field and the orientation of the molecule, which 
inspires its name: scalar coupling (also called indirect spin-spin coupling, or J-coupling). 

B. The NMR Hamiltonian 

As mentioned above, in liquid state NMR the intermolecular spin interactions arc sup- 
pressed. This causes the molecules to be efi^ectively isolated from each other, and therefore 
a description of the spin dynamics of an ensemble of molecules is well approximated by the 
spin dynamics of a single molecule. If we consider a molecule containing N spin-1/2 nuclei 
in one of the molecules, then the natural Hamiltonian of this system in a large homogeneous 
magnetic field Bq pointing in the z direction is given by 

1 ^ 

i=l i<j 

where i^f = cof^/^Ti = 7i|-Bo| is the Larmor frequency of the i^^ nucleus with gyromagnetic 
ratio 7j, Jij is the coupling strength between nucleus i and j and cr* is the z Pauli matrix of 
the i^^ spin. 

The first term in the Hamiltonian describes the precession of the spins due to their 
coupling to the external magnetic field, while the second term describes the J-coupling 
between pairs of nuclei. This Hamiltonian corresponds to the weak coupling limit, where we 
assume that the chemical shifts between coupled spins are much greater than their respective 
couplings, i.e. {uj^ — » Jjj/2. If this approximation is not valid, we need to use the 
full coupling Hamiltonian a*.a;^ + o'^o'^ + c^o"^ in place of ct*(T^. The exact values of the 
Hamiltonian parameters are determined by fitting experimental data. 

C. Single-spin control 

For quantum information processing, we need to be able to perform arbitrary manipu- 
lations of a single spin, which is equivalent to arbitrary rotations about any axis. As an 

example, consider the application of a magnetic field Bi perpendicular to the z axis which 
oscillates at the nuclear spin's Larmor frequency: 

Si = I Si I (cos(a;^-'^ t)x + sm{uj'^^ t) y) (2) 
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where = 2tiv'^^ is the angular frequency of the field. In the rotating frame of the nucleus 
(i.e. the frame rotating at the same frequency as the spin), Bi will appear as a constant 
field pointing along its rotating x axis. The spin will start to precess about this axis. 
Rotation about any axis in the xy-plane is possible by adjusting the phase of the Bi field, 
e.g. uj^t — > uj^t + (f), which will create a rotation around the axis making an angle (f) with the 
X axis. In the laboratory, such a rotating field can be applied by sending a radio-frequency 
(RF) pulse of a particular duration and phase to a conducting coil surrounding the sample, 
calculated according to the rotating wave approximation (see [4] for more details). 

To better understand this phenomenon from the viewpoint of quantum mechanics, con- 
sider the rotating frame picture: suppose the spin is in the state \%lj{t)), and define the state 
in the rotating frame of the pulse with angular frequency uj'^^ as 

= i?,(-cu"^t)e~^^"-|7/;(0)) 

= es'"^V*e"7?^"'"|^(0)) (3) 
= |'?/'(0)), for a single spin with a; '^•^ = a;^. (4) 

If we apply a time derivative to equation [3} it can be shown that the state in the rotating 
frame evolves according to the Schrodinger equation with the new Hamiltonian 

rf 

Hr = R.{-Uj'h)nnatRzW'h) - ^^a,. (5) 

When an RF pulse with phase is applied to the spin, the laboratory frame Hamiltonian 
is: 

n = —(T, + — (cos {uj'H + 0)cT^ + sin {uj'H + (j))ay) (6) 

where a;™* = 7r7i|Si|. In the rotating frame this becomes 

Hr = ^(cu^ -cu'^O^^. + ^^""*(cos0a^ + sin0(T3,). (7) 

Thus, if the RF pulse is at the same frequency as the spin, the spin will see a constant field 
in the xy plane, and will precess about it. The rotation angle 6 is determine by the interval 
T during which the RF field is applied, according to ^ = cj^V. 
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D. Adding a second spin 

It is also possible to independently control two spins with different Larmor frequencies. 
Applying an RF pulse at the frequency of the first spin, the rotating frame Hamiltonian is 
given by 

Unat = ^^rV^ + + \{u^, - u{)al + lj,2alal (8) 

where we have set = for simplicity. While the first spin undergoes a rotation around 
the X axis, the second spin experiences a field with an additional non-zero z component. 
This is called the off resonance effect. If we consider the case where — >> cj""* then 
the second spin rotation around the x axis will average to zero during the time the first 
spin has completed its rotation. Typically, o;"^"* is smaller than 1 MHz, so this condition is 
automatically satisfied if the two nuclei belong to different species. If the spins are of the 
same species, this condition can also be satisfied if a very low amplitude pulse is used due 
to the small nutation frequency. In this case, one drawback is that the pulse will necessarily 
take much longer to achieve the same angle of rotation, and if the two spins have a significant 
coupling constant J12 coupling effects might introduce significant errors and therefore limit 
our control. 

Fortunately, there exist well known techniques to address different nuclei of the same 
species with high precision. The most common technique is to control the spins using 
shaped pulses. The frequency response to the pulse will depend on the pulse shape (Fourier 
theorem) and so by applying the pulse with a time varying power we can control the power 
spectrum of the pulse. For example, if a Gaussian shaped pulse is applied at frequency a;''-^, 
then only spins within a Gaussian distribution of frequencies around uj"^^ will respond to this 
RF field. Therefore, if the height and the length of the Gaussian pulse is carefully chosen, 
one spin can be "addressed", causing negligible effects to others. This technique permits 
control of spin pairs with smaller chemical shift differences in shorter periods of time, hence 
allowing stronger coupling. The length of a Gaussian pulse is proportional to the inverse of 
the chemical shift between the spins. Therefore, in the limiting case of small chemical shift 
differences and large J-coupling values, control of the qubits is more difficult. 

For most liquid state experiments on a few spins, where chemical shifts are comparatively 
large and J-couplings are small, the use of Gaussian pulses is sufficient to achieve very 
high precision spin rotations. The situation becomes more complicated when there are 
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more homonuclear spins (implying smaller chemical shift differences on average) , or stronger 
coupling like in solid state or liquid crystalenvironments. It is still possible to overcome 
these drawbacks by considering more complicated pulse shapes and phase modulation. For 



example, in section |I G[ we will describe how it is possible to find shaped pulses that can 
implement any desired evolution by simulating the full quantum dynamics. 



E. The controUed-NOT operation 

In the previous subsection we discussed a method used to independently control different 
spins. In order to perform quantum computing, we need to achieve universal control and 
hence be able to have spins interact with each other. A two-qubit gate that is useful for 
quantum information processing is the controUed-NOT, which acts as 

|00) ^ |00) , |10) 111) 

|01)^|01) , |11)-|10). (9) 

The operation must flip the target qubit (second bit) if and only if the flrst qubit is in the 
state |1). In NMR, |0) and |1) are associated with the state of the spin pointing up, | t) or 
pointing down | J.) respectively. If we look at the Hamiltonian in equation ([T]), and consider 
its effect on spin 2 depending on whether the state of spin 1 is up or down, we obtain an 
effective Hamiltonian for the second spin: 

H] = ]^{uj!^ + TiJi2)(yl (10) 

nl = \{u;^-'nJr,)al (11) 

Therefore, spin 2 will rotate slower or faster depending on the state of spin 1. If the coupling 
evolves for a time r = 2j^, we obtain a controlled-Z rotation of | degrees, which can be 
transformed into a controlled-NOT by a few single spin pulses applied before and after (see 
flgure[l]for the complete sequence to implement a controlled-NOT). 

In the previous sections, we demonstrated that it is possible to implement any rotation 
around an arbitrary axis in the xy-plane, as well as perform a controlled-NOT gate with 
two spins. These two conditions are sufficient to perform universal quantum computing. 
For practical purposes, it is convenient to use only ^ and vr RF pulses as they are easier 
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FIG. 1: Implementation of a controUed-NOT gate in liquid state NMR. The left circuit is the 
quantum circuit representation of a controUed-NOT gate with control qubit qi and target qubit 
q2. On the right is the NMR implementation of such a gate by combining single qubit rotations 
and the natural two spin interaction of the system. The single qubit rotation properties are given 
in the rectangles, e.g. X : 90 is a short notation for Rx{^) = e~*4'^^. Notice that the size of the 
rectangles are not to scale for liquid state NMR, where the RF pulses are usually much shorter 
than the time to implement a | J-coupling 



to calibrate. Arbitrary angle rotations can still be applied because 2;-axis rotations can be 
obtained by changing the definition of the rotating frame, which is equivalent to changing 
the phase of subsequent pulses. It can be verified that 



(12) 



where the vector n = cosOx — sm6y. Therefore, since z rotations commute with the 
internal Hamiltonian of the system, we may commute all the z rotations to the end of the 
pulse sequence and compensate for any remaining z rotation during the post-processing 
of the data. Moreover, the overall z-rotation takes no time and is far more precise when 
using this procedure because the RF phase has a higher accuracy than the RF amplitude in 
modern NMR spectrometers. 



F. Refocusing and control techniques 



In an NMR system, spins constantly couple to each other, and we must "turn off" these 
couplings on demand to implement generic quantum gates. For example, consider a three 
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spin system in which we wish to implement a controlled-NOT between the first and second 
qubits. As mentioned above, a | J-couphng between spins 1 and 2 is needed, which is 
accomphshed by allowing the system to evolve under the natural Hamiltonian for a time 
J—. During this time spin 3 will also couple to spins 1 and 2, giving an unwanted 



r 



2Jl2 



evolution. However, if we apply a vr pulse on the third spin half-way through the free 
evolution (at time |), this spin will effectively decouple from the other two spins and, upon 
an extra tt pulse at the end of the evolution (at time r), it will be returned to its initial 
state (this pulse sequence is shown in Figure [2]). Considering only the interaqction term of 
the intramolecular component of the Hamiltonian, we can write the evolution of the system 
as: 

= (^(j3,)e~*T(-^12'^z'^^+'^13'^z'^z+-^23f7zCTf ^_^^3'jg-i2f (Ji20-1ct2^ 
— g-«i(-'l20-lCT2-Ji30-lo-|-J23(T2CTf)g-i5{Ji20-lcr2 + Ji3CT>| + J23(T2crf) 

where -R^(vr) is the operator of a 7r-pulse about the x-axis on spin 3. This is called a 
refocusing scheme and can be readily generalized to any number of coupled spins, i.e. a vr 
pulse on spin i will effectively decouple it from all the other spins. This scheme can also be 
efficiently generalized to decouple m spins from each other and from the N — m remaining 
spins. In practice, the situation is more complex. For example, one must keep track of the 
Zeeman evolution of all the spins (which is called phase tracking). This evolution can be 
taken into account by changing the phase of subsequent pulses according to the relation 



given in equation 12 



For systems involving up to a few spins, pulse phases and decoupling sequences are derived 
by hand, but for molecules containing greater number of spins, these calculations become 
tedious and computer assisted techniques are used [5]. Efficient classical algorithms can 
be implemented that optimize pulse sequences with respect to phase and residual coupling 
errors. 

A major source of pulse errors are off resonance and coupling errors. It is possible to 
estimate and compensate for these errors by evaluating the first-order coupling and phase 
errors generated by a pulse. This is done by assuming that the real pulse can be decomposed 
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FIG. 2: In NMR, the coupling between spins is always active. It is possible to refocus 2-qubit 
interactions using special pulse sequences. An example is given above, where halfway through a 
period r, a vr pulse on one of the nuclei is implemented that reverses the direction of its spin. Note 
that when we have a coupling of the form (f^^^^ff^^^^ the effect of the pulse can be thought of as 
reversing the sign of the coupling, and thus allows to cancel the interaction that occurred during 
the first r/2 period. This pulse sequence effectively decouples the third qubit from the system, 
while leaving the coupling between qi and q2 unchanged. This situation is mathematically treated 



in equation (13) 



as the ideal pulse preceded and followed by phase and coupling errors (see figure [3]). Since 
the error terms all commute with each other, they can be estimated using pairwise spin 
simulations, which requires reasonable computational resources, i.e. is efficient as we scale 
the number of nuclei. With small J-couplings and short pulses it is reasonable to expect 
error rates below a fraction of a percent for each pulse. 

Once the errors generated by each pulse are known they can be taken into account and 
corrected for by optimizing the durations of the free evolution periods and the timing of the 
refocusing pulses. Such an algorithm can also perform phase tracking and modify the pulse 
phases accordingly. Very high gate fidelities have been demonstrated using this efficient 
pairwise simulation technique. 
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FIG. 3: A selective pulse designed to implement a single-qubit rotation in an A^-spin system will, 
in general, also affect the other spins. This can be studied in small systems by simulating the full 
quantum dynamics to obtain the unitary Ugim- The unwanted evolution of the non-target qubits 
is represented by the broken line on the left figure. If the pulse is carefully designed so that its 
implementation is very close to the ideal desired unitary Uid, the error can be efficiently estimated 
to first order by phase errors (represented by Oi and (pi) and coupling errors occurring before and 
after the pulse. 

G. Optimal control for strongly coupled spins 

In some cases, spins are so strongly coupled that the approximation of crlcri couplings 
used above breaks down. In those cases another technique can be used: strongly modulating 
pulses designed using numerical optimal control techniques [HI [7]. For systems containing 
about less than ten qubits, we can find extremely high fidelity and robust control by applying 
optimal control principles. Just as classical optimal control theory can tell how to best steer 
a rocket, quantum optimal control gives the tools to best steer a quantum system to a 
desired unitary gate. Quantum optimal control has been used for some time in the context 
of driving chemical reactions with shaped laser pulses [S]. There, the goal is to maximize 
the transfer from a known initial state to a known final state. In the context of quantum 
computing, we do not necessarily know what the input state will be, and so we must find 
unitary gates which will work correctly for any input state. 

The Hamiltonian at any point in time can be written down as 



(14) 
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where, Tinat is the time-independent, natural (or drift) Hamiltonian, and T-Lrf{t) represents 
the controllable time-dependent RF fields which we can use to drive the evolution of the 
system. The task is to find the sequence of control fields that will produce the correct 
unitary evolution. For reasonably complicated systems analytical solutions are beyond reach; 
however, using a numerical search we can find control sequences with fidelities as high as 
0.999999. The control fields are discretized at a suitable rate (see figure |4]) and a random 
guess for the control sequence is chosen. We can simulate the evolution of the system under 
this sequence and obtain the simulated unitary Usim- We can then compare this to our goal 
unitary Ugoai through a fitness function. 



which is equivalent to the state overlap fidelity averaged over all input states [6]. We can 
then use any optimization routine to search for the highest $. 

The first step in this direction for NMR quantum information processing was the scheme 
of Fortunato et al. [6J who, by limiting the form of the control fields to a small number 
of constant amplitude, phase and frequency periods, were able to find high fidelity control 
sequences through a simplex search. Many more control periods can be considered and 
the numerical search substantially sped up through the use of standard optimal control 
techniques to obtain information about the gradient of the fitness function at each point. 
This is the GRadient Ascent Pulse Engineering or GRAPE algorithm introduced by Khaneja 
et al. [7J. From the simulation information we can calculate approximate gradients of the 
fitness function with respect to the control amplitude at each timestep (see Figure |4]). With 
this gradient information we can update the control fields by moving along the steepest 
ascent direction and then repeat the procedure. Convergence to a global maximum is of 
course not guaranteed as this hill climbing algorithm is a local search. However, empirically 
we find we can achieve sufficiently high fidelities from such local maxima. Convergence of 
the algorithm can be further improved by using non-linear conjugate gradients [H]. 

The control sequences drive the system through a complicated and non-intuitive path and 
small errors in our modeling of the system's Hamiltonian might lead to a low fidelity pulse. 
Fortunately, the control sequences can also be made robust against static inhomogeneities or 
uncertainties such as field inhomogenties and amphtude miscalibration of pulses. Robustness 
to both these effects for a particular pulse is shown in figure [6]d. 




(15) 
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FIG. 4: An example section of a GRAPE pulse. The bars show the control amplitude which is 
constant for each time step. The arrows show the derivatives from the GRAPE algorithm at each 
point which tell us how to update the pulse for the next iteration (adapted from [7]). 

Recent work on a 12-qubit liquid state system that compared this method with the more 
traditional approach outlined in previous sections can be found in |TUJ. Not surprisingly, 
the control using strongly modulating pulses was more precise but at an exponential cost in 
searching for the pulse sequences. 

The scheme described above can be generalized if we can break the molecule into strongly 
coupled subsystems which are themselves weakly coupled, and thus it is possible to perform 
the decomposition described in previous sections between subsystems using pairwise simula- 
tion of subsystems. Gates to implement rotations within the subsystem will however require 
numerical optimization methods such as GRAPE. The details of this hybrid method are 
presently under investigation. 

The ability to control liquid state nuclear spin systems has allowed implementation of a 
variety of benchmarking experiments and algorithms on small qubit registers. For example, 
we have developed sufficient control to implement quantum error correcting protocols pTl 
\TI\ [131 [3 [IS]) the simulation of quantum systems ^1^1 [THl [IHl 120] and the benchmark 
12 qubit quantum processor [TD]. An extensive list of these and other experiments can be 
found at http:/ /arxiv.org/. 
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II. SOLID STATE NMR 

Solid state NMR presents a very different spin environment in comparison to liquid state 
NMR. Dipolar couplings do not average to zero in solids and in fact tend to dominate 
the NMR spectra, usually through significant broadening of the resonance peaks. As an 
illustration, consider an array of identical spins on a crystalline lattice: any one nuclear spin 
couples most strongly to its nearest neighbors, then to next-nearest neighbors, and so on 
(the dipolar coupling is long range, falling off as where r is the distance between spins). 
At low spin polarizations, each coupling induces a splitting of the spin's resonance energy 
into two components of slightly higher and lower energies. Averaging over the ensemble 
of spins, a smoothly broadened resonance is observed, whose width is comparable to the 
size of the nearest-neighbor coupling. These couplings can be quite large compared to the 
J-couplings observed in liquids; for example, a proton pair separated by 2A has a maximal 
dipolar coupling frequency of 15 kHz, whereas proton-proton J-couplings are usually 200 
Hz or less. 

The solid state spin environment would at first appear to be less than ideal for quantum 
information processing, as the strong spin-spin interactions would lead to rapid effective loss 
of coherent quantum information. However, this situation can be remedied in two ways: (1) 
the evolution due to time-independent interactions is, in principle, completely refocusable by 
appropriate RF pulse sequences and (2) the spatial dilution of spins in the sample can be used 



to minimize spin-spin interactions. For the three-qubit processor described in section II A 
based on the malonic acid single crystal, we take advantage of both spin dilution and NMR 
refocusing techniques to maintain good qubit coherence. Additionally, spin-lattice relaxation 
times Ti can be very long (minutes to hours) in solids, particularly for spin-1/2 nuclei in 
crystalline solids at low temperatures and with low impurity or defect densities. We generally 
have Ti > T2 » T2, where is the dephasing time in the absence of refocusing, and T2 
is the "intrinsic" coherence time which is usually determined in practice for solid state 
systems by the quality of the refocusing controls. With perfect control, we should obtain 
T2 = 2Ti, which would allow for many thousands of quantum gates to be performed within 
the characteristic coherence time scale. Of course, in practice, control errors and other 
experimental imperfections greatly limit the effective spin coherence time. 

Besides the potential for much longer coherence times, solid state NMR also offers the 
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prospect of greatly increasing spin polarization by cooling the sample to low temperatures. 
Nuclear Zeeman energies at attainable magnetic fields are, however, on the order of 1 /zeV or 
less, and such small energies require a temperature less than 10 mK to reach full polarization 
by thermal equilibration. Fortunately, nuclei may also be polarized dynamically at much 
more accessible temperatures ~1.5K, by driving polarization exchange between nuclear spins 
and a small concentration of unpaired electrons that have been introduced into the sample 
[21] . Since the magnetic moment of the electron is about 660 times larger than that of the 
proton, electron spins can be fully polarized in thermal equilibrium at such temperatures. 

More generally, nuclear and electron spins incorporated into a solid lattice are thought 
to be promising elements for building scalable quantum information processors [221 
[21]. Nanofabrication technologies are continually growing in precision and sophistication, 
allowing ever more controllable nanoscale devices. Nuclear spins can provide stable qubit 
storage, and the transfer of information to electron spins can allow for fast operations and 
coupling to optical photons can provide initialization, readout and long-range entanglement 
[231 123 [2S1 [2Z] • Controlling the nuclear and electron spins with high fidelity using magnetic 
resonance should play an important role in these schemes. 



A. Three qubit solid state quantum information processor 

The following example of a small quantum information processor is based on malonic 
acid, C'iHiOi [29]. The sample is a single crystal (grown from aqueous solution) consisting 
of a mixture of ^^C^H^O^ and ^^C^H^Oi (labeled) molecules; ^^C is spin-1/2 whereas ^^C 
is spinless. The fraction of labeled molecules is kept to < 10% of the total, so that the 
^^Ca molecules can function as an ensemble of 3-qubit registers with relatively weak spin 
interactions between registers. It is important to note that all molecules contain 4 spin-1/2 

atoms, so that generally, a decouphng RF sequence must be applied to the spins to 
isolate the ^^C qubits from the spin system. Figure ?? shows a ^^C NMR spectrum (with 

decoupling) obtained at an external field magnitude and orientation at which the three 
carbons are well separated in Larmor frequency (i.e. separately addressable) and the carbon- 
carbon intramolecular dipolar couplings are relatively large, allowing for fast coupling gates. 
The spin Hamiltonian of the 3-qubit system is determined directly by fitting a simulated 
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FIG. 5: (right) Fitted spectrum of a malonic acid single crystal (3.5% ^^C3-ff404 concentration) 
at a particular orientation with respect to the external field, (left) Table showing the Hamiltonian 
parameters extracted from the spectral fitting: chemical shifts are along the diagonal, dipolar 
couplings in the lower off-diagonal (green), and J-couplings in the upper off-diagonal (yellow). 
The dipolar parameters involving the methylene proton spins Hmi,2 have been estimated indirectly 
by using the ^^C parameters to precisely determine the crystal orientation, and knowledge of the 
crystal structure from neutron scattering data [28] to then calculate the couplings. 

spectrum to the measured one. The height differences between the absorption peaks are due 
to the "strong coupling" effect, in which non-diagonal terms in the coupling Hamiltonian 
(such as duio'lo'l + (^y^y)) have frequencies di2 comparable to the difference in chemical 
shifts of the coupled spins {/S.uj^{a], — a1)). 

When di2 ~ Au;^ as in the present solid state system, the a\ eigenstates (i.e. the 
computational basis states |000), |001), etc.), are no longer the energy eigenstates and the 
spin dynamics become more complex. This poses a challenge for universal quantum control: 



the scalable method described in section |I F| for generating quantum gate pulse sequences 
in liquid state NMR is valid only when such mixing is negligible. However, numerically 
optimized strongly modulating pulses can succeed in providing sufficient control. 

This control can be made robust with respect to sources of ensemble incoherence. For 
example, the linewidths ((^2*)^^) of the resonance peaks visible in Figure ?? are about two 
orders of magnitude larger than those seen in a liquid environment. This incoherence is 
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FIG. 6: (a) An example of a GRAPE pulse designed to implement a controlled-NOT-NOT on the 
three qubit malonic acid processor, (b) The pulse is designed to be robust to inhomogeneities in 
RF amplitude and external magnetic field; the fidelity of the pulse is plotted versus both types of 
inhomogeneity. The blue curve shows the fidelity of the gate averaged over all input states while 
the red curve shows the fidelity for the worst-case input state. Note that the robustness of the gate 
fidelity to static field deviation is equivalent to self-refocusing of T| dephasing by the pulse. 



mainly due to a combination of residual dipolar couplings between ^^CsH^O^ molecules and 
inhomogeneity of the external magnetic field over the sample [32] ■ Strongly modulating 
pulses can be made sufficiently robust to these effects with minimal added computational 
overhead, particularly when using the GRAPE algorithm (see figure [6]). 



B. Algorithmic Cooling 

As a computational demonstration of the malonic acid quantum processor, we now con- 
sider the experimental implementation of a heat-bath algorithmic cooling (HBAC) protocol 
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FIG. 7: (a) Schematic quantum circuit diagram of the HBAC protocol implemented in malonic 
acid. Time flows from left to right. The last step corresponds to a three-bit compression unitary 
operator, and is here decomposed into two controUed-NOT-NOT gates sandwiching a Toffoli gate 
|30j . (b) Experimental results, in terms of measured qubit polarizations at each step. The full 
bars indicate the ideal values of the protocol, and the shaded bands give the experimental values 
(uncertainties are given by the widths of these bands). 



|31j . The aim of the experiment is to amplify the polarization on one of the three qubits, 
using both qubit-qubit operations and a controlled coupling to a surrounding "heat-bath" 
[321 E3] ■ Note that HBAC is a classical algorithm, since each step results in a classical state 
of the system, however we use quantum gates in its implementation. Strongly modulating 
pulses obtained using GRAPE are used to perform gate operations on the qubits. The HBAC 
algorithm operates on the three malonic acid qubits as well as the large surrounding spin 
system which acts as a spin-bath with a large effective heat capacity. Figure [7^ shows a 
schematic quantum circuit diagram of the first six steps of the HBAC protocol implemented 
in the malonic acid system. The swap operations between Cm and Hmx (labeled 'refresh' in 
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figure [Tp,) can be carried out either using a dual ^^C /^H multiple-pulse refocusing sequence 
[311 [3l] or by using Hartmann-Hahn cross-polarization [35] with short contact time |40j . 

The multi-qubit operations are carried out using strongly modulating pulses. During 
these operations, a strong continuous wave decoupling field is applied to the spins, 
serving both to decouple them from the carbons and to "lock" their collective magnetization 
along a particular direction in the rotating frame (the latter is called 'spin-locking' and is 
used to prevent the magnetization from dephasing). This continuous wave field permits 
the protons to interact with each other via the — ^ H dipolar couplings. Such couplings 
diffuse spin polarization [36] through the Hamiltonian terms of the form cr^o"^ + a\(j\ [JT]. 
Therefore, the Hmi polarization of the spins restore to the equilibrium value of the 
bath after sufficient interaction time. The timescale for such equilibration in this system is 
~ 50 X T;{H) ^ 5ms « Ti{H). 

Experiments were performed at room temperature in the regime of low nuclear spin polar- 
ization, and the aim was to demonstrate quantum control similar to what has been achieved 
using liquid state NMR quantum information processors. The measured polarizations on 
each of the qubit spins at each step of the algorithm are shown in Figure [7j The full bars 
indicate ideal values and the shaded bands give the experimental values. A rough indication 
of the fidelity of the experimental protocol can be gauged by the ratio of the final polariza- 
tion of the target qubit C2 to its ideal value, F ^ 1.39/1.5 = 92.7%. This yields an average 
error (here loss of polarization) per step of about 1.5%. Given that the inhomogeneous line- 
broadening effects (i.e. the natural incoherence of the system) in solid state NMR is two 
orders of magnitude larger than that in liquid state, this is an impressive degree of control. 
Furthermore, these results indicate that such sequences applied to single quantum systems 
would yield much higher fidelities. We are presently implementing a multiple-round HBAC 
algorithm with the goal of moving beyond the 1.5 signal amplification benchmark that uses 
only unitary transformation. 

III. CONTROL OF COUPLED ELECTRON-NUCLEAR SYSTEMS 

An attractive direction for spin-based quantum information processing is to couple the 
nuclear spin degree of freedom with that of the electron spin. The larger magnetic moment 



20 







\ 






ViV. 


— ITi) 








Vl3 


V23' 


— ^1 U) 


3 1 




>|2co„| 









Vis 



1 V / 



(d) 

2 



FIG. 8: (a) The ENDOR model for universal quantum control of an electron-nuclear hyperfine 
coupled spin pair. All states can be reached using three control parameters: two electron tran- 
sitions at microwave frequencies (blue arrows) and one nuclear transition at radio-frequency (red 
arrow), (b) If the system has an anisotropic hyperfine interaction, universal control can be achieved 
using only one microwave control field. The mixing of states connected by "forbidden" transitions 
(dashed arrows) occurs under the hyperfine interaction, allowing pathways to all states. Mod- 
ulated sequences of the single control field can be found using the GRAPE method to generate 
arbitrary unitary gates, (c) Stable radical of malonic acid, produced by x-ray irradiation. One of 
the methylene protons is removed by irradiation, leaving behind an unpaired vr-electron spin-1/2. 
The maximum hyperfine couplings to the methylene ^^C and are 220 and 90 MHz, respectively, 
(d) A representative level connectivity diagram for a le-2n system in an orientation for which the 
anisotropic hyperfine terms are large. Universal quantum control is achievable by implementing 
appropriate modulated sequences on the microwave transition labeled ui^. A ^^C labeled nucleus 
on the methylene carbon of malonic acid would exhibit a similar level structure. 
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of the electron leads to a much faster timescale for spin manipulation, and the electron's 
charge degree of freedom can serve as a pathway for initialization and readout of spin qubits. 
Early progress in this direction was made by Mehring et al. [371 EH], who showed coherent 
control of an ensemble of electron-proton spin pairs in irradiated malonic acid using electron- 
nuclear double resonance (ENDOR) techniques. This control scheme is sketched in FigurejS^, 
where universal quantum control is achieved in principle by pulse sequences involving two 
microwave and one RF field. 

Moreover, Mehring introduced the important idea that, in hyperfine systems with one 
electron and > 1 distinguishable nuclear spins, quantum gates operating on nuclear states 
could be sped up significantly using the hyperfine interaction with the electron to mediate 
indirect nuclear-nuclear interactions [24j. This was dubbed the S-bus concept since the 
electron spin operator is often denoted as S. The S-bus concept is suggestive of a picture 
in which isolated clusters of 1 electron + nuclear spins form qubit sub-registers, where a 
handful of qubits can be manipulated and stored with high fidehty. 

Initialization of the electron spins to the ground state |0) occurs thermally at modestly 
low temperatures and typical external magnetic fields. The nuclear spin qubits can be ini- 
tialized one by one by alternating polarization swaps with the electron spin and waiting 
for subsequent electron thermalization (Tf''^'^ << T""'^). There are various possibilities for 
coupling sub-registers: the electron-electron exchange coupling modulated by electrostatic 
gates [22], direct dipolar coupling between electrons, or long-range coupling using optical 
techniques [231 [2S1 12SI I2Z]- A recent proposal outlines a method to further improve the com- 
putational power of the S-bus concept. This proposed method uses the natural anisotropy 
of the hyperfine interaction to "mix" spin sublevels (much as the strong coupling effect dis- 
cussed above mixes nuclear spin computational basis states). This eliminates the need for 
direct RF control of the nuclear spin transitions. This method is shown schematically in 
Figure |8)d. The effective spin Hamiltonian of the electron-nuclear pair has the form: 

Hen = UJ^S, + UJ^h + A,S,h + A,S,h, (16) 

where A^. {Az) is the anisotropic (isoptropic) hyperfine coupling, S and I are the electron 
and nuclear spin operators, and uo^^ are the Larmor frequencies. If we arrange the system 
such that Ax ~ uj^., there will be efficient mixing of the sub- levels indicated by the dashed 
arrows in figure [8]d, due to the hyperfine term A^SzIx- 
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It can be shown mathematically that controlled driving of a single microwave transition 
(indicated as the 1-3 transition in figure [sjo) is enough to generate any unitary operator 
when given sufficient bandwidth, duration and complexity of the control field. This remark- 
able result allows for much faster operations by replacing the low energy nuclear excitation 
field with the higher energy hyperfine coupling. It also simplifies the experimental setup by 
requiring only one control frequency. Moreover, we can apply optimal control pulse-shaping 
methods such as GRAPE to generate universal quantum gates. Experimental work is cur- 
rently underway with the aim of demonstrating these techniques using hyperfine coupled 
spin systems such as the stable malonic acid radical (see figure Isb-d). 



IV. CONCLUSION 



Of the many proposed and current realizations of quantum information processing de- 
vices, NMR stands out as the present leader in terms of control fidelities and the manipu- 
lation of the largest number of quantum bits. This has allowed for the implementation of a 
variety of benchmarks and algorithms that bring theoretical concepts to the laboratory. Its 
success is built on many years of RF engineering invested in NMR technology for a variety 
of purposes that have recently been adapted for quantum information processing. Combin- 
ing such control with the long coherence times of nuclear spins makes NMR promising for 
quantum computation. 

One of the present difficulties with using only nuclear spins is the ability to initialize them 
into pure states. Therefore the next generation of magnetic resonance based devices will 
make use of electron spins both to initialize nuclear spin qubits and to mediate interactions 
between them, yielding much faster gate operations. We expect this research to continue to 
produce exciting results and to merge with other quantum technologies in the coming years. 
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